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ABSTRACT 

The chemical properties of dwarf spheroidals in the local group are shown to be inconsistent with star 
formation being truncated after the reionization epoch (z~8). Enhanced levels of [Ba/Y] in stars 
in dwarf spheroidals like Sculptor indicate strong s-process production from low-mass stars whose 
lifetimes are comparable with the duration of the pre-reionization epoch. The chemical evolution 
of Sculptor is followed using a model with SNe II and SNe la feedback and mass- and metallicity- 
dependent nucleosynthetic yields for elements from H to Pb. We are unable to reproduce the Ba/Y 
ratio unless stars formed over an interval long enough for the low-mass stars to pollute the interstellar 
medium with s-elements. This robust result challenges the suggestion that most of the local group 
dwarf spheroidals are fossils of reionization and supports the case for large initial dark matter halos. 

Subject headings: galaxies: abundance ratios — galaxies: dwarf — galaxies: Local Group — galaxies: 
evolution — stars: abundances — stars: AGB — nucleosynthesis 



1. INTRODUCTION 

Dwarf galaxies in the Local Group provide a nearby 
laboratory for testing dark matter (DM) structure for- 
mation scenarios and indirectly probing the epoch of 
reionization. Many properties of Milky Way's dwarf 
spheroidal (dSph) satellites have been well studied, yet 
their role in the cosmological context is still uncertain. 
Cosmological simulations have long been plagued by the 
"missing satellites" problem - the overprediction by up to 
two orders of magnitude of the abundance of DM subha- 
los with mass < 10 8 M in systems like our Local Group 
(Klypin et al. 1999). The two most popular explana- 
tions for the observed deficit of low-mass satellites are 
that: 1) only massive (> 1O 9 M ) satellites were able to 
form stars, with the spheroidal morphology and small 
present-day stellar velocities of dSphs arising through 
tidal interactions (e.g. Kravtsov et al. 2004; Kazantzidis 
et al. 2004); or 2) low mass satellites (< 10 8 M Q ) could 
form stars, but only prior to the reionization epoch, af- 
ter which UV heating photoevaporates the gas and cur- 
tails star formation (Barkana & Loeb 1999; Bullock et al. 
2000; Ricotti & Gnedin 2005; Gnedin & Kravtsov 2006). 
These two solutions have very different, and testable, 
consequences for the dSph population. The first scenario 
implies that Galactic dSphs descended from far more 
massive DM halos than one infers assuming that mass 
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follows light (i.e. 10 7 ~ 8 M Q , Mateo 1998), while the sec- 
ond scenario implies that many dSphs are very ancient 
(z > 8) and associated with less massive DM halos. 

Recent work from Ricotti & Gnedin (2005) favors the 
second scenario, in which most dSphs are "fossils" of 
reionization. They suggest that internal ionizing sources 
and, to a lesser extent, UV background radiation pre- 
vented the gas from cooling and collapsing to form 
stars after the reionization era. In this case, almost all 
star formation (SF) activity transpired in an interval of 
< 600 Myr 1 . This relatively brief window of opportunity 
for SF should leave observable signatures on the stellar 
chemical properties. A commonly used "clock" of SF 
is a/Fe versus metallicity (Matteucci & Recchi 2001). 
This is because a-elements like O, Mg and Si are gener- 
ated primarily in short-lived massive stars culminating as 
Type II supernovae (SNe II), whereas a substantial frac- 
tion of Galactic Fe is attributed to Type la SNe (SNe la), 
which occur on longer timescales (Timmes, Woosley & 
Weaver 1995). The first generation of SNe II pollute 
their host galaxy with supersolar a/Fe gas, while a de- 
layed contribution from SNe la subsequently leads to a 
decrease in a/Fe with increasing Fe/H. Thus, a system 
with low a/Fe at low metallicity is interpreted as having 
undergone slow and protracted SF, whereas stars born in 
a short intense burst should exhibit higher a/Fe. Unfor- 
tunately a/Fe does not uniquely trace SF. Dwarf galaxies 

1 here we have assumed ACDM cosmology with present-day 
Mm = 0.34, uj\ = 0.66, Hubble parameter = 70 km s _1 Mpc - 1 , 
and a 13 Gyr age for the universe 
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are especially vulnerable to feedback effects that may fos- 
ter the preferential loss of elements from massive SNe II, 
thus complicating the interpretation of the a/Fe diag- 
nostic. 

In this paper we employ the abundance pattern of 
neutron-capture elements, particularly those associated 
with the slow neutron-capture process (known as the 
s-process), as a better chemical test of SF duration in 
dSphs. Venn et al. (2004) suggested that the overabun- 
dance of [Ba/Y] in most dSphs was consistent with the 
production of s-elements from lower-mass stars. We con- 
firm that the curious behavior of the relative abundance 
of neutron-capture elements Y, Ba and Eu in dSphs re- 
flects a strong contribution from low-mass stars, which 
are the chief source of the s-elements (Busso et al. 1999). 
Owing to the correspondingly long lifetimes of these low- 
mass stars, one may use chemical evolution simulations 
to place limits on the minimum duration of star forma- 
tion allowed by the observations. In Section|2we give the 
details of our chemical evolution model of the represen- 
tative dwarf galaxy Sculptor. Section [3] presents results 
from models with different star formation histories, in 
order to assess whether dSphs are survivors or fossils of 
the epoch of reionization and Section^Jcontains our main 
conclusions. 

2. THE MODEL 

We have modeled the chemical evolution of the dwarf 
spheroidal Sculptor on the assumption that its evolution 
and SF history is typical of the types of dSph galaxies 
that Ricotti & Gnedin (2005) identified as likely fossils 
of reionization. Moreover, its proximity makes Sculptor 
one of the best studied dSphs with a large set of chemical 
abundance constraints (e.g. Shetrone et al. 2003; Geisler 
et al. 2005). This section summarizes the key features of 
the model, which is described in greater detail in Fenner 
et al. (2006, in prep), where it is applied to ten local 
group dwarf spheroidals and irregulars. 

2.1. Equations 

The chemical enrichment history of Sculptor was sim- 
ulated using a modified version of the chemical evolution 
code described in Fenner & Gibson (2003) and Fenner et 
al. (2003, 2004, 2005). Major changes made to the code 
include: allowing for SNe-driven outflows; employing the 
star formation history from Dolphin et al. (2005); and 
adding s- and r-process yields. The evolution of the gas 
phase abundance pattern reflects the cumulative history 
of the dynamic processes of star formation, gas inflows 
and outflows, stellar evolution, and nucleosynthesis. The 
set of equations governing these processes were numer- 
ically solved by defining ai{t) as the mass surface den- 
sity of species i at time t, and assuming that the rate of 
change of <7j(t) is given by: 



— = Ei ,LIMs(t) + Ei SNIl(t) + Ei sNIa{t) 

at 

-Wi JSM (t) - W it SNIl(t) ~ Wi,SNIa(t) 

+ ^<ri{t)infall-Xi[t)l/j(t), (1) 

where the first three terms, Ei^iMs(t), Ei^Niiit) and 
Ei t sNia(t), denote the mass surface density of species 



i ejected at time t by low- and intermediate-mass stars 
(LEVIS), SNe II and SNe la, respectively. The follow- 
ing term, WijsM(t), denotes the mass surface density 
of species i that escapes from the neutral galactic in- 
terstellar (ISM) due to being carried along with SNe- 
driven winds. The next two terms, Wi t sNii(t) and 
Wi^Nia{t), represent mass of species i released by SNe II 
and SNe la, respectively, that is expelled via winds. The 
term -^Ui(t)i n f a u gives the infall rate of i, which is deter- 
mined by calculating the amount of gas required to en- 
sure the SFR satisfies the Kennicutt law (see Section l2~2"|l . 
The chemical composition of the infalling gas (and the 
gas at the first timestep) was taken to be primordial. 
The final term, Xi(t) ip(t), gives the depletion of species 
i due to incorporation into newly formed stars, where 
Xi(t) is the mass fraction of i in the ISM at time t and 
tf)(t) is the SFR. The ejection rates in EqnQJ specifying 
the production of each element from the different types 
of star, are given by: 
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f m " v 4>{m) 
i,SNii(t)= / ip(t-T m )Yi(m, Z(t-T m )) dm, 

Jm bu m 

(4) 

where Yi(m, Z(t — r m )) is the stellar yield of i (in mass 
units) from a star of mass m, main-sequence lifetime r m , 
and metallicity Z(t — r m ). The initial mass function, 
4>(m), has upper and lower mass limits of m up =60MQ 
and m/ oto =0.08 Mq, respectively. The models presented 
here employ the Kroupa, Tout & Gilmore (1993) three- 
component IMF. The fraction of LIMS in binary systems 
culminating in a SNe la is given by the parameter A. 
The SNe la progenitor binary mass is denoted by m& 
and is the sum of the primary and secondary masses, 
mi +TO2. The binary mass has upper and lower limits of 
mi m =16M and 771^=3 M Q , respectively. The function 
f(n) determines the distribution for the mass fraction of 
the secondary (li = m 2 /mfc). The reader is referred to 
Matteucci & Greggio (1986) for more detailed discussion 
and definitions of the terms in the above equations. 

Most previous studies (e.g. Lanfranchi et al. 2003; 
2004; 2006; Robertson et al. 2005) have assumed dif- 
ferential galactic winds in proportion to the SFR. These 
have assumed Wi(t) = Wiip(t), where Wi is a free pa- 
rameter that is arbitrarily set for each element in or- 
der to reproduce observations, and is then kept constant 
as a function of time. Such a scheme does not allow 
for behavior such as increased escape of Fe relative to 
O after a burst of SF when the frequency of SNe la 
may compare with, or exceed, that of SNe II. Instead, 
we self-consistently calculate the amount of SNe II and 
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SNe la ejecta lost at each timestep as a function of the 
appropriate SNe rates. Cold interstellar gas also escapes 
from the galaxy by imposing a mass-loading factor that 
is consistent with observations. The form of the equa- 
tions governing the loss of species i in galactic winds can 
be written: 



W hS Nia(t) = Ei. S Nia{t) xmin(0.9, £sNIa xR SN Ia) (5) 

mtot 

W iiS Nii(t) = E itSN ii{t) x mm(0.9, x R S nii) (6) 

mtot 



W i>ISM (t) = Xi(t) X ML X (W gaSt SNIa(t) + W gaSlS NIl(t)) 

(7) 

where m to t is the total galaxy mass in units of 10 6 M Q 
(taken from Mateo 1998). The Type la and II SNe rates 
are given by RsNia and Rsnii, while the correspond- 
ing wind efficiency factors are Esnio, and Ssnii- The 
mass of element i in the ISM that gets entrained in the 
SNe-driven winds and driven out of the galaxy is set by 
the mass-loading factor, ML (= 15, as discussed below). 
The sum of W gas ,SNla(t) and W gaSt sNii{t) gives the to- 
tal amount of stellar ejecta lost in winds at time t, and 
is equivalent to the sum of Wi t sNia{t) and Wi t sNii{t) 
over all i. The escape fraction of SNe II and SNe la 
ejecta increases linearly with their respective SNe rates, 
but the overall efficiency is inversely proportional to the 
depth of the galactic potential well, as measured by the 
total galaxy mass, mtot- We do not allow more than 90% 
of the SNe ejecta to escape at each timestep, based on 
the prediction by Hensler et al. (2004) that up to 25% 
of SNe material gets mixed into the local ISM on short 
timcscales rather than being carried away by an expand- 
ing superbubble. We note that the main conclusions of 
this paper are robust to the implementation of feedback. 

Martin et al. (2002) estimated that the diffuse halo 
gas of the dwarf starburst NGC 1569 has a/Fe ~ 2- 
4 times solar. The metallicity of the hot halo gas was 
measured to be roughly solar, implying a mass-loading 
factor of -9. In order to obtain [O/Fe] ~ +0.3 to +0.6 
in the galactic winds, we assume that SNe la and SNe II 
have the same explosion energy of 10 51 erg and we set 
the feedback efficiency for SNe la to be five times higher 
than that of SNe II. This might seem counter intuitive, 
given that SNe II occur in associations, which should en- 
hance their ability to heat the ISM and create bubbles 
and chimneys facilitating the escape of metals. How- 
ever, Recchi et al. (2004, 2006) argue that SNe II are 
less efficient at expelling their ejecta, because much of 
their energy goes toward heating the cold dense molecu- 
lar clouds from which they were born. The characteristic 
lifetimes of SNe la range from tens to thousands of Myr, 
by which time the SNe la progenitors will have traveled 
to regions of lower density, where less energy is required 
to widely disperse ejecta. 

We assume that turbulent motions at the boundary 
between the SNe-driven winds and the ambient ISM mix 
significant quantities of cold gas into the outflows. Our 
adopted mass-loading factor of 15 is similar to the value 
inferred for NGC 1569 and the value of 10 calculated 
by Silk (2003). The mass- loading factor (i.e. the mass 
of the ISM carried away with the winds relative to the 
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Fig. 1. — Star formation rate in units of Mq per thousand years 
as a function of time for models A and B. 



mass of stellar ejecta in the wind) is the key parameter 
determining the metallicity of the galactic outflows. 

Observations of superwinds in starburst galaxies by 
Hcckman et al. (2000) reveal outflow rates (OFR) com- 
parable with the SFR. Theoretical calculations support 
this approximate equivalence between the OFR and SFR 
(Silk 2003; Shu, Mo & Mao 2005). Silk (2003) esti- 
mated that OFR ~ SFR, in the case of large porosity, 
and OFR < SFR, when porosity is low. In accordance, 
our overall wind efficiency was chosen to produce outflow 
rates roughly half the SFR. We also assumed that galac- 
tic wind efficiency is inversely proportional to the total 
galaxy mass, where galaxy mass was taken from Mateo 
(1998). 

Given a velocity dispersion of about 2 km/s for stars 
in OB associations (Tian et al. 1996) and lifetimes 
> 60 Myr for low- and intermediate-mass stars (Schallcr 
et al. 1992), then most LIMS will have have drifted out- 
side the ~ 50-100 pc radii of the stellar association in 
which they were born by the time they begin to pollute 
the ISM with their winds. Owing to their weak stellar 
winds and their likely migration from OB bubbles, we 
mix the newly released ejecta from LIMS directly into 
the ISM at each timestep. 

2.2. Star Formation History (SFH) 

We tested two different star formation histories: 1) 
model A adopts the SF history inferred from Sculptor's 
color-magnitude diagram (Dolphin et al. 2005) in which 
SF extends for over several Gyrs; 2) model B compresses 
the SF within a period of 600 Myr. Fig. □ compares 
these two SF histories. The gas accretion rate is often 
treated as an input into chemical evolution models, how- 
ever we work backwards to infer the infall rates from the 
SFHs by inverting the Kennicutt (1998) SF law. Like 
Lanfranchi et al. (2003, 2004, 2006), we find that low 
SF efficiencies are better able to reproduce the chemi- 
cal abundances in local group dwarf galaxies. In partic- 
ular, we obtain good fits in our standard model using 
ip(t) = 0.05 x ff(t) L4 M pc~ 2 Gyr -1 , where ip is the star 
formation rate (SFR) and a is the gas surface density. 
The difference between the o~(t) implied by the tp(t) and 
the value predicted within the chemical evolution mod- 
els, provides the infall rate at each timestep. 
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2.3. Stellar Yields 
LOW AND INTERMEDIATE MASS STARS (LIMS) 

For the production of isotopes up to the Fe-peak from 
stars less massive than 8 , we incorporated yields from 
the stellar evolution and nucleosynthesis code described 
in Karakas & Lattanzio (2003), supplemented with un- 
published yields for the metallicity Z=0.0001. The origin 
and interpretation of s-process elements like La, Ba and 
Y in the local group is of considerable import to this 
study. The Karakas et al. models do not calculate the 
s-process, so for elements above the Fe-peak, we employ 
the nucleosynthetic yields of the Torino group (Gallino 
et al. in prep). These mass and metallicity dependent 
yields have previously been used in the Galactic chemical 
evolution models of Travaglio et al. (1999, 2001, 2004). 
In a similar fashion to Travaglio et al., we have averaged 
the s-process yields over a range of 13 C pocket choices. 

The slow neutron-capture process in asymptotic giant 
branch (AGB) stars is thought to be driven primarily 
by the C(a, n) 16 reaction, which provides the flux of 
neutrons that can be captured by elements from the Fe- 
peak upwards (Busso et al. 1999). The source of the 13 C 
probably comes from CN-cycling in the H envelope lead- 
ing to proton capture by 12 C and the development of a 
pocket of 13 C at the top of the intershell region between 
the H-shell and He-shell. This pocket arises through the 
diffusion of a small amount of protons from the H enve- 
lope at the epoch of the third dredge-up episodes. Proton 
mixing is induced by the close contact between the con- 
vective H-rich envelope and the radiative He-rich and C- 
rich He intershell, this forming a tiny "proton pocket". 
These protons are subsequently captured by the abun- 
dant 12 C in the He intershell, creating a 13 C pocket of 
primary nature. Afterwards, when the pocket is heated 
up to about 0.9 xl0 8 K, 13 C is consumed by a captures 
that release neutrons and give rise to a very intense neu- 
tron exposure. The formation of the 13 C pocket is not 
simulated within the nucleosynthesis codes, but rather is 
treated parametrically and constrained by observations 
(Busso et al. 1999). 

TYPE IA SUPERNOVAE (SNE IA ) 

We adopted a recalculation of the Thielemann, 
Nomoto, & Yokoi (1986) W7 model by Iwamoto et al. 
(1999) to estimate the yields from SNe la. It was as- 
sumed that 4% of binary systems involving intermediate 
and low mass stars result in SNe la, since this fraction 
provides a good fit to the solar neighborhood (e.g. Alibes, 
Labay, & Canal 2001; Fenner & Gibson 2003). 

MASSIVE STARS 

For stars more massive than 8-10 M that end their 
lives in violent supernova explosions, we implemented 
the yields from Woosley & Weaver (1995) for elements 
up to Zn (using same technique as Fenner et al. 2003, 
2004 and halving the Fe yield, as suggested by Timmes, 
Woosley, & Weaver 1995). Beyond Zn, the yields of r- 
process elements like Eu are difficult to calculate within 
nucleosynthesis models and it is common practice to as- 
sume primary production from SNe II and empirically 
deduce them (e.g. Argast et al. 2004) or estimate them 
by calculating the difference between the solar abun- 
dance and the pure s-process contribution predicted by 
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Fig. 2. — Predicted abundances for the Sculptor dSph chemical 
evolution model, incorporating SN feedback and galactic winds, 
plotted against stellar observations from Shetrone et al. 2003 and 
Geisler et al. 2005 (diamonds). The shaded region indicates the 
relative frequency of stars from model A as a function of [X/Y] 
and [Fe/H], where [X/Y] is indicated in the top lefthand corner of 
each panel. Pairs of contour lines are plotted at 0.25 and 0.75 of 
the maximum frequency. The predicted stellar abundances have 
been convolved with a Gaussian of dispersion 0.1 dex to mimic 
observational uncertainties. Results from model B, in which SF is 
shut-off by reionization after 0.6 Gyr, are shown by dashed contour 
lines. Dotted contour lines correspond to model C, which is the 
same as model A but without SNe feedback. 



Galactic chemical evolution models (e.g. Travaglio et 
al. 1999, 2001; 2004). Although the stellar mass range 
responsible for r-production is still under debate (e.g. 
Sumiyoshi et al. 2001; Wanajo et al. 2003; Tsujimoto et 
al. 2000), our models link the r-process site to all SNe II 
with mass < 40 M Q . We tested the case in which only 
SNe II with mass < 15 M Q produced r-elements, finding 
that the conclusions of this paper are not very sensitive 
to the adopted mass range. Metallicity-dependent main- 
sequence lifetimes calculated by Schaller et al. (1992) 
have been employed. 

3. RESULTS 

Fig. El shows the predicted abundances of O, Mg, Si, 
Ca, Sc and Ti relative to Fe from the Sculptor model 
using: i) the standard SFH inferred from observations 
(model A: solid contour lines and shaded region); and 
ii) the suppressed SF model (model B: dashed contour 
lines). To illustrate the role of feedback and preferential 
loss of SNe ejecta, we also show results from an identical 
simulation to model A but without galactic winds (model 
C: dotted contour lines). Red giant stars observed in 
Sculptor by Shetrone et al. (2003) and Geisler et al. 
(2005) (diamonds) have a mean [Fe/H] ~ —1.57, which is 
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Fig. 3. — Same as Fig. 121 but for neutron-capture elements. The 
shaded region indicates the relative frequency of stars from model 
A as a function of [X/Y] and [Fc/H], where [X/Y] is indicated in the 
top lefthand corner of each panel. Pairs of contour lines are plotted 
at 0.25 and 0.75 of the maximum frequency. The predicted stellar 
abundances have been convolved with a Gaussian of dispersion 
0.1 dex to mimic observational uncertainties. Results from model 
B, in which SF is shut-off by reionization after 0.6 Gyr, are shown 
by dashed contour lines. Dotted contour lines correspond to model 
C. 



well matched by our predicted mean of [Fe/H] ~ — 1.65. 
We also produce a metallicity spread in agreement with 
the data. Sculptor stars exhibit a /Fe ratios 2 slightly 
above solar at the lowest metallicities, before declining 
to subsolar values at higher [Fe/H] (>— 1.5). Note that 
the subsolar [a/Fe] values in Sculptor stars contrasts dra- 
matically with the supersolar ratios characterizing MW 
stars at the same metallicity. This is thought to be 
common in dwarf galaxies whose low gas surface den- 
sities lead to inefficient SF and consequently long metal 
enrichment timescales. However, we find that subsolar 
a/Fe can arise through the preferential loss of SN prod- 
ucts via galactic winds (compare models A and C). Both 
the standard and suppressed SF models satisfactorily re- 
produce the overall a-element trends. The standard SF 
model has slightly lower (by ~ 0.1 dex) O, Si, Ca, Sc and 
Ti with respect to Fe than in the suppressed SF model, 
but the influence of galactic winds makes the difference 
smaller than for a closed box model. The prevention of 
SNe feedback (dotted contour lines) leads to an excess of 
O, Si and Ca, in conflict with the data. Fig. [21 demon- 
strates that [a/Fe] is not a clean probe of SF history in 
galaxies that have experienced metal-enriched outflows. 

2 where a-elements include O, Mg, Si, Ca and Ti. Although we 
also plot Sc in this figure, it is classed as an iron-peak element. 



Indeed, [a/Fe] is more sensitive in our Sculptor models 
to SNe feedback than to the duration of SF. 

Conversely, Fig.[3]suggests that the neutron-capture el- 
ements provide a more useful discriminant of the different 
star formation models. The neutron-capture elements 
show curious patterns in Sculptor stars (e.g. Venn et al. 
2004). The light s-process element Y tends to be subso- 
lar relative to Fe, and lower than in Milky Way stars of 
comparable metallicity. Conversely, the heavy s-process 
element Ba is slightly higher than in Galactic field stars. 
In our models, the Sculptor stars with enhanced Ba/Y 
ratios were born from gas that had been heavily pol- 
luted by low-mass stars. The high [Ba/Y] ratio can only 
be matched by our model incorporating both SNe-driven 
winds and SF extending over several Gyr. This is be- 
cause the origin of gas with [Ba/Y] > +0.5 dex is low- 
mass metal-poor stars on the asymptotic giant branch, 
whose weak (~ 15-30 km s _1 ) winds are easily retained 
by the potential well. 

For low-mass AGB stars to dominate s-element abun- 
dances in the ISM there needs to be substantial loss of 
SNe II ejecta and SF timescales exceeding the lifetimes 
of low-mass stars. The strong metallicity-dependence of 
s-process yields from AGBs gives rise to the sharp up- 
turn in Ba and La abundance at [Fe/H] ~— 1.6. This 
leads to some stars having Ba-enhancements as high as 
observed, although the fraction of simulated stars with 
supersolar Ba relative to Fe, Y and Eu is underestimated, 
even in model A. This may inform us about the size and 
nature of the 13 C pocket in metal-poor AGB stars. AGB 
stars produce heavier s-elements like La and Ba more 
efficiently than lighter species like Y at low metallici- 
ties. This is because the neutron flux per seed nuclei 
varies roughly inversely with metallicity. Since the slow 
neutron-capture process proceeds further along the peri- 
odic table with decreasing initial stellar metallicity, Pb 
production from AGBs peaks at lower [Fe/H], followed 
by a Ba and La peak at higher [Fe/H]. The production 
of Y from AGBs reaches peak efficiency at higher [Fe/H] 
than is typical for dSphs. A further prediction from our 
model is that dwarf spheroidal stars should have super- 
solar [Pb/Ba] and [Pb/Fe] ratios. Although Pb requires 
very high signal-to-noise and has not yet been measured 
in stars outside the Milky Way, sufficient exposure time 
may yield future detections in red giants in nearby dSphs. 
Such measurements would be a very useful test of this 
model and may provide insight into the nature of the 
s-process in low-mass stars. 

The relative importance of the s- and r-process is 
usually measured by the [Ba/Eu] ratio, where the pure 
r-process ratio is [Ba/Eu] ~ —0.8 dex and the pure s- 
process ratio is ~+ldex (Arlandini et al. 1999). Fig- EI 
shows about an order of magnitude rise in Ba/Eu go- 
ing from [Fe/H] = — 2 to —1. Sculptor stars tend towards 
[Ba/Eu] = +0.5 at [Fe/H] = -1, which is 0.8dex above 
the corresponding MW value. Only Model A reproduces 
this trend, owing to the retention of lower-mass stellar 
ejecta and the preferential loss of the ejecta from SNe II, 
which are major suppliers of Eu. The good agreement 
between Model A predictions and the observed [Ba/Fe], 
[La/Fe] and [Ba/Eu], lends further support to the case 
for extended SF in dSphs. We place more weight on the 
Ba/Y diagnostic, however, since there is no obvious alter- 
native to low-mass stars that can match the high values 
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Fig. 4. — [Ba/Y] in the ejecta from AGB stars versus initial 
stellar mass for the three metallicities indicated (for a choice of 
13 C pocket mass of ~ 4 X 10 — 6 Mq) Corresponding stellar lifetimes 
from Schaller et al 1992 (for Z=0.001) are indicated along the top 
axis. 



measured. 

Interestingly, the massive Galactic globular cluster oj 
Centauri exhibits similar chemical properties, with the 
lowest metallicity stars ([Fe/H~ —1.8) having a predom- 
inantly r-process abundance pattern and those with the 
highest metallicity ([Fe/H] ~ — 0.8) having dramatically 
enhancements in s-process elements (Smith et al. 2000). 
Smith et al. attributed these observations to the greater 
retention of low-mass stellar ejecta - because of their 
weak winds - compared with energetic SNe II winds. 
They infer that star formation must have continued for at 
least several Gyr, making u Centauri highly anomalous 
compared with other globular clusters, whose individual 
stellar populations formed contemporaneously. Indeed, 
it has been suggested (e.g. Majewski et al. 2000) that u> 
Centauri is not a globular cluster, but the remains of an 
accreted dwarf spheroidal. 

To further illustrate the time delay associated with s- 
process enrichment, Fig. 0] displays the ratio of [Ba/Y] 
in the ejecta of AGB stars as a function mass for the 
three metallicities [Z/H]=— 1.9, —1.6 and —1.2 for the 
standard choice of 15 C pocket mass (see Travaglio et 
al. 1999, 2001, 2004 for further details regarding the 
13 C pocket). The corresponding stellar lifetimes (from 
Schaller et al. 1992, for a metallicity of Z=0.001) 
are presented along the top axis. Sculptor stars have 
a mean <[Fc/H]> ~-1.6dex and < [Ba/Y] > 0.5 dex. 
From Fig. 01 it can be seen that at metallicities typi- 
cal of dSphs, only AGB stars with mass <3M Q and 
lifetimes >300Myr have yields with [Ba/Y] > 0.5 dex. 
Stars whose yields have [Ba/Y] equal to the highest val- 
ues measured in Sculptor stars have lifetimes > 1 Gyr. 
If these stars polluted the ISM with Ba-rich gas from 
which subsequent generations of stars formed, then stars 
must have been forming for several Gyr or more. We 
note that with a factor of 1.5 reduction in the 13 C 
pocket, slightly more massive AGB stars can yield 
[Ba/Y] > 0.5 dex on shorter timescales (mass < 4 M Q and 
lifetimes >150Myr). Indeed, smaller 13 C pockets may 
help resolve the underproduction of [Ba/Y] and [Ba/Fe] 



seen in Fig. [3J After weighting by the IMF, however, 
lower-mass stars dominate the production of s-elements 
and still require extended SF in order to pollute the ISM 
to the observed levels. 

We have argued that the suppression of SF by reion- 
ization after ~ 600 Myr would not grant low-mass stars 
enough time to pollute the ISM with s-elements to the 
extent demanded by the observations. However, there is 
still uncertainty regarding the precise redshift at which 
the universe was reionized. To test the sensitivity of our 
results to the assumed redshift of reionization, we ran an 
identical model with z re i on ~ 6, which is roughly the low- 
est value permitted by observations (e.g. Benson et al. 
2005, Fan et al. 2005). This model generously allowed 
SF to continue for 0.9 Gyr, yet the mean [Ba/Y] and 
[Ba/Eu] increased by less than 0.15 dex and still violated 
the empirical constraints. 

Grebel & Gallagher (2004) found little evidence for 
reionization suppressing local group galaxy evolution, 
based on: i) metallicity spreads; ii) low [a/Fe] ratios; and 
iii) stellar ages. However, a range in metallicites along 
with low [a/Fe] can both arise from a short star forma- 
tion episode (see dashed contour lines in Fig. 0). The 
clearest sign of ongoing SF would be an age spread, how- 
ever, age-dating old stellar populations is complicated by 
the age-metallicity degeneracy and is subject to uncer- 
tainties > 1 Gyr (e.g. Hernandez et al. 2000). Using 
s-element enrichment from long-lived low-mass stars as 
an alternative chronometer, we arrive at the same con- 
clusion as Grebel & Gallagher (2004). 

If early reionization was not responsible for shutting off 
SF in dSphs like Sculptor, what ultimately caused star 
formation to cease after a few Gyrs? In our models, SNe 
feedback can remove the bulk of the metals from Sculp- 
tor (and the other dSphs, as shown in Fenner et al. 2006, 
in prep) but only a small fraction of the ISM gets car- 
ried away by the SNe-driven winds. Consequently, our 
Sculptor model has a final gas fraction of ~ 70%, whereas 
most dSphs are observed to have gas fractions less than a 
percent. Lower gas fractions could be obtained by invok- 
ing unreasonably strong galactic winds, however a more 
likely mechanism for gas removal is via tidal- and ram- 
pressure stripping when the orbits of the dSphs bring 
them into close proximity with the Milky Way halo (e.g. 
Marcolini et al. 2003, 2004, 2006; Mayer et al. 2005). 

4. DISCUSSION AND CONCLUSIONS 

We find that the abundance of neutron-capture ele- 
ments in dwarf spheroidal stars is indicative of a strong 
contribution from long-lived low-mass stars. While this 
study presented models of Sculptor, our conclusions ap- 
ply equally to other old gas-poor dSphs like Draco, Ursa 
Minor and Sextans, since they also exhibit enhanced 
Ba/Y and have similar star formation histories to Sculp- 
tor (Dolphin et al. 2005). The observed supersolar Ba/Y 
ratios imply that stars formed over an interval of at least 
several Gyr to allow time for metal-poor AGB stars to 
enrich the ISM with s-elements. These results cast doubt 
on recent suggestions that local group dSphs Draco, Ursa 
Minor and Sextans are "fossils" of reionization (Ricotti 
& Gnedin 2005; Gnedin & Kravtsov 2006). Because the 
Ba/Y ratio is very high in all dSphs for which it has 
been measured, we propose that all dSphs were actively 
forming stars beyond the epoch of reionization. This is 
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indirect evidence for their initial dark matter halos being 
larger than the ~ 10 7 M Q values obtained by assuming 
that mass follows light (Mateo 1998). 

The case for large mass-to-light ratios and extended 
DM halos in dSphs continues to strengthen. Mashchcnko 
et al. (2005) estimated virial masses of ~ 10 9 M for 
Draco, Sculptor and Carina, based on N-body simula- 
tions designed to fit the observed stellar density and ve- 
locity profiles. Dchncn & King (2006) indirectly infer 
a DM halo mass > 10 9 M Q for Sculptor, from its large 
abundance of X-ray binaries, whose high expected ve- 
locities should have facilitated their escape if the po- 
tential well was smaller. Moreover, there appears to 
be an absence of tidal heating at large radii in Draco, 
Ursa Minor, Sextans and Sculptor, which implies large 
and extended DM halos (Read et al. 2006; Coleman et 
al. 2005). The chemical properties of the representa- 



tive dwarf spheroidal Sculptor are well matched by our 
simulations incorporating SNe-driven feedback and em- 
pirical star formation histories. We contend that the 
neutron-capture abundance pattern reflects enrichment 
from metal-poor low-mass stars, which is inconsistent 
with a cessation in star formation after the reionization 
epoch and instead supports dSphs having large initial 
DM halos. 
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